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CO Abstract 

We present a quantitative analysis of the effect of rough hydrophobic surfaces on 
£Nj viscous newtonian flows. We use a model introduced by Ybert and coauthors in |17j . in 

(— i which the rough surface is replaced by a flat plane with alternating small areas of slip 

and no-slip. We investigate the averaged slip generated at the boundary, depending on 
the ratio between these areas. This problem reduces to the homogenization of a non-local 
C — ■ system, involving the Dirichlet to Neumann map of the Stokes operator, in a domain 

with small holes. Pondering on the works [21 [3J, we compute accurate scaling laws of the 
averaged slip for various types of roughness (riblets, patches). Numerical computations 
p I complete and confirm the analysis. 

< 

^ 1 Introduction 

o3 

With the development of microfluidics, drag reduction for low Reynolds number flows, notably 
at solid walls, has become a stimulating issue. Therefore, the interaction between a fluid and 
a solid boundary has been investigated thoroughly, both at the experimental and theoretical 
levels. A special attention has been paid to the detection of slip, for various types of flows 
and solid walls. We refer to |13| for a review. 

As a result of this activity, the idea that rough boundaries could generate a substantial slip 
has spread out. This idea has developped on the basis of both experimental and theoretical 
works, notably on wall laws. We remind that in the context of roughness effects, a wall law 
CO is an effective boundary condition imposed at a smoothened boundary, reflecting the overall 

impact of the real rough boundary. In particular, if one describes the rough boundary through 
an oscillation of small amplitude and wavelength e, one can show rigorously that a no-slip 
condition at the rough boundary can be replaced by a wall law of Navier type, with slip length 
5— i of order e. We refer for instance to [U [T2l [5] for more precise statements. 

a 

However, these seemingly favorable results must be considered with care. For instance, 
at the experimental level, one must ensure that the slip is not measured too much away 
from the boundary. Also, as regards the theoretical works on wall laws, the position of the 
artificial boundary at which the law is prescribed is crucial. In particular, all forementioned 
works consider artifical boundaries that are at the top of the roughness. As a result, the flow 
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rate in the smoothened domain does not equal the averaged flow rate in the rough domain, 
making comparisons inaccurate. This is the issue of the so-called apparent slip, noticed in []. 
In fact, in the case of rough wetting surfaces (endowed with a no-slip condition), one can even 
show the following: if one puts the artificial boundary in a way that the flow rates are the 
same, then the flat boundary is optimal with respect to drag minimization. We refer to [6] for 
detailed statements and proofs. Hence, the possibility of decreasing drag through roughness 
is not so clear, especially for rough wetting surfaces. 

Still, in the recent years, promising results have been obtained concerning a class of rough 
hydrophobic surfaces, see for instance [16J. Indeed, by the combination of the chemical and 
geometrical properties of these surfaces, the hollows of the roughness get filled with gas. 
Hence, the viscous fluid above does not penetrate: it slips above the hollows, and only sticks 
at the bumps, reaching the so-called Cassie or fakir state. 

The aim of this paper is to study the slip generated by such configurations, both in a 
rigorous and quantitative manner. We focus on a model proposed in article [17], in which the 
rough boundary is replaced by a flat plane, divided in small periodic cells (say of side e -C 1). 
Each cell is divided in two zones: 

• A no-slip zone, corresponding to a plane projection of the sticky part of the roughness 
(bumps). 

• A slip-zone, corresponding to a plane projection of the slippery part. 

Using homogenization techniques, we derive an effective boundary condition as e goes to zero, 
depending on the size a € of the no-slip zones. We provide in this way scaling laws for the slip 
coefficients, for various configurations (patches, riblets). Such laws are in global agreement 
with the formal computations led in |17] . One shows notably that the riblet configuration is 
less effective than patches one (see Remark^. All our theoretical results are grounded by 
numerical computations at the end of the paper. 



2 Main results 

Let us first present the model under study. We consider a three-dimensional Stokes flow 
between two infinite plates: 

-An + Vp = f in 

(2.1) 

div u = in f2, 

where VL = T 2 x (0,1) and T 2 = M 2 /Z 2 . We denote by x = (x\, X2, £3) = (^,£3) the space 
variable. The function / £ L 2 {Q) is a given source term. On the upper surface X3 = 1, we 
enforce a "no-slip" boundary condition 

"1*3=1 = 0. (2.2) 

On the lower surface, we assume that u satisfies alternately "perfect slip" and "no slip" bound- 
ary conditions, corresponding respectively to the hollows and bumps of the rough hydrophobic 
surface. More precisely, let e > and 

5 e := [0,6) 2 ~ (M/(eZ)) 2 , 
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the elementary square of side e . For simplicity, we shall assume all along that e _1 is an 
integer. Let T € be a Lipschitz subdomain of S e , modeling an elementary no-slip zone. Details 
about T € will be given right below. From this elementary no-slip zone, we define a global one 
inside [0, l) 2 ~T 2 : 



[J (ek + T e 



fee [lo,. 



Finally, the boundary condition at X3 = is 

-"31x3=0 = 0, d3U h \ X3=0 = on (T e ) c x {0}, 1*^=0 = 0onfx {0}. 



(2.3) 



It is easily proved that J27I^-([2T2|)-([2T3|) has a unique solution (u e ,p e ) G H 1 ^) x L 2 (tt)/R. 

This article is devoted to the asymptotic analysis of (u e ,p e ), as e — > 0. We will distinguish 
between two types of no-slip pattern T e : 



Patches: we assume that 



2' T 



+ a e T, 



(2.4) 



where (f , |) is the center of the square S e , and where the domain T is relatively compact 
in (—1/2, 1/2) 2 , and contains a disk of radius a > 0, centered in the origin (see Figure 
[TJ. The parameter a e is a positive number such that a e < e. In this case, the no-slip 
zone is a union of periodically distributed patches. 



Riblets: we assume that 



T e := (eT)x( 2+ aJ). 



(2.5) 



where I C (— |, |) is an open interval (see Figure ^b. In this case, the no-slip zone is a 
union of stripes, invariant in the xi-direction. Of course, invariance in the X2 direction 
could have been considered as well. Note that later on, addressing the case of riblets, 
we shall focus on two particular cases: 

— / = ey. riblets parallel to the flow; 

— / = &2- riblets perpendicular to the flow. 



The issue is to derive a wall law for the system (2.1)-(2.2)-(2.3), that is, to replace the mixed 



boundary condition (2.3) at X3 = by a condition which does not depend on e. We will show 



that u e behaves asymptotically like the solution u in H 1 of d2T|-((2T2l), endowed either with 
a Navier boundary condition 



«3 = at X3 = 0, d-^Uh = Muh at X3 = 0, 



or with a Dirichlet boundary condition 



u 



|x 3 =0 



0. 



(2.6) 



(2.7) 



In ( 2.6 ), M is a 2 x 2 non-negative matrix, whose eigenvalues have the dimension of the inverse 
of a length. If M = Aid, the number A -1 is called the "slip length". In the general case, 
the inverse of the eigenvalues provide the slip lengths in the directions of the eigenvectors. 



We shall denote um the solution of (2.1 )-(2.2)-(2.6). We will write uq in the special case 
M = (0 o)- Eventually, we shall denote Uoo the solution of (p7T])-(f2T2|)-(p?f|) . 



With the previous notation, we can state our first result: 
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Figure 1: Patch configuration. For every k = (k±, k 2 ) £ [|0, . . . , e" 1 ]] 2 , the intersection of the 
no-slip zone T e with the cell [eki, e(k\ + 1)) x [ek 2 , e{k 2 + 1)) is defined by ek + T e = ek + a e T. 
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Figure 2: Riblet configuration. For k = (ki,k 2 ), the intersection of the no-slip zone T e with 
the cell [efci, e(fci + 1)) x [ek 2 , e(k 2 + 1)) is defined by ek + T £ = ek + (eT) x (§ + a e 7). 
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Theorem 1. (Asymptotic behavior for patches) 

Assume that T e := (§,§) + a e T. Let u e G H 1 ^ 2 x (0, 1)) be the solution of ([21]), §2?2\ , 



(2.3). One must distinguish between three cases: 

1. Sub-critical case: if a t <C e 2 , then u € — uq in H 1 ^ 2 x (0, 1)); 

2. Super- critical case: if a t 3> e 2 , then u € — -Soo in ff 1 (T 2 x (0, 1)); 

3. Critical case: there exists a symmetric, positive definite matrix Mq such that if a e /e 2 — > 
C > 0, then u € u Co m in H 1 ^ 2 x (0, 1)). 

A similar result holds for riblets. Let us merely state the theorem in the critical case: 

Theorem 2. (Asymptotic behaviour for riblets) 

Assume that T € := (eT) x (a e I). Suppose that lim e _*.Q — eln(a e ) = Co > 0, and furthermore 
that f does not depend on x\ . 

Then, u e — 1 u u rib , where 

M ri6 =(f I). (2.8) 
\ Co J 

Additionally, when f = e\ or f = e2, the limit system can be simplified: 

• if f = e\ (riblets parallel to the main flow), then UM rib ,2 = UM rib ,3 = and u~M rib ,i 
satisfies 

dzu Mrib ,i = ^UMribA at x 3 = 0. 
Hence, the slip length is Cq/tt; 

• if f = e2 (riblets perpendicular to the main flow), then u~M rib .i = and UM rib ,2 satisfies 

2vr_ 

03U Mrib ,2 = yr u M rib ,2 at X 3 = 0. 

Hence, the slip length is Co/(2tt). 

Remark 3. Notice that theorems^ and\^ indicate that for a given area of no-slip, riblets are 
less efficient than patches in terms of slip optimization. Indeed, consider riblets T e with area 
ea € corresponding to the critical case in Theorem [s| lim^o — eln(a e ) = Cq > 0. Patches with 
the same area have a diameter of order ^ea t : they have a sub-critical behavior (with regards 
to Theorem^ since y / eal <C e 2 . So, u e converges towards u$. In other words, for the same 
area, patches achieve perfect slip on the lower surface, while riblets only have finite slip. In a 
similar fashion, it can be proved that riblets whose area a 2 is critical for patches - i. e. such 
that a e ~ e 2 , see Theorem^ have a super- critical behavior (with regards to Theorem^) and 
converge towards Uoo-' the solution satisfies a no-slip condition at the limit. 

Remark 4. Our results are consistent with those of [T 7 ^ : indeed, in the case of patches, it 
is shown heuristically there that the slip length is proportional to e 2 /a e : in other words, if 
a e < e 2 , perfect slip is achieved, if a e ^$> e 2 , a no-slip condition is retrieved in the limit, and 
in the critical case, the slip length is positive and finite. 

Also, explicit calculations (see WM) recalled in show that the slip length for riblets 
is equal to — e/7rln(a e /e) for riblets parallel to the flow, and to —e/(2ir)\n(a e /e) for riblets 
perpendicular to the flow. Once again, this is consistent with Theorem [1| 
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Remark 5. Theorems 1 and 2 do not support the idea that rough hydrophobic surfaces can 
generate a substantial slip. Indeed, to obtain an effective slip law, the surface fraction of 
no-slip has to be very small. Back to wall roughness, this would correspond to narrow peaks 
separated by (comparatively) large hollows. It seems far from the roughness characteristics 
used experimentally to obtain a hydrophobic Cassie state. 

The proofs of theorems 1 and 2, that rely strongly on the papers [H E] by Allaire, are 
given in Section [3] and [4] respectively. We then present in Section [5] numerical simulations 
that confirm the asymptotic results, and clarify the influence of the shape of patches on the 
slip length, i.e. on the eigenvalues of the matrix M. 



3 Asymptotic study of "patch" designs 

This section is devoted to the proof of Theorem [l] Let (u e ,p € ) G H 1 ^) x L 2 (Q)/R be 



the solution of (2.1), (2.2), (2.3). By classical arguments, the sequence (u e ,p e ) is uniformly 
bounded in x L 2 (0)/K, and consequently, there exists a couple (u,p) G x 

L 2 (tt)/R such that 

u e weakly in p e —^p weakly in L 2 (Q)/M.. 



Using the weak formulation of Eqs. (2.1 ) and the continuity of the trace operator, one obtains 



easily that the weak- limit (u,p) satisfies Eqs. (2.1) and boundary condition (2.2) on x% = 1. On 



X3 = 0, the boundary condition satisfied by the vertical component is preserved in the limit, 
and we obtain u^ X3= q = 0. To describe the boundary condition satisfied by the horizontal 
components (ui,U2) on X3 = 0, we need to distinguish between the so-called super-critical, 
critical and sub-critical cases. 

Notation. For every k = (/ci,/^) £ [ 1 , e 1 1 ] 2 , we denote the elementary squares, cubes 
and half-cubes as follows: 

S e k := ek + S\ P£ := 5| x (-|, |), P^ + := n R%. (3.1) 

We shall use that notation throughout the paper. 

Super-critical case: a e 3> e 2 . The proof in the super-critical case relies on a quantitative 
Poincare inequality: we claim that there exist eo > and a positive function 77(e) such that 
77(e) — > as e — > 0, and such that 



f K| 2 < 77(e) f \Vu e \ 2 V0<e<e . (3.2) 
jt 2 x|o| J a 



/T 2 x{0} 

We provide a proof of this inequality in the Appendix. 

Since u e is bounded in ff 1 (f2), we immediately infer that T^a^o vanishes in L 2 (T 2 ). Thus 
(u,p) is a solution of the Stokes system with homogeneous Dirichlet boundary conditions at 
X3 = and X3 = 1, i.e. u = Uoq. 

Critical and sub-critical cases: a e < e 2 . We follow here the strategy of articles [21 13] 
by Allaire. These articles deal with the homogenization of the Stokes equations across a 
network of balls, with a Dirichlet condition at the surface of the balls. Notably, in section 
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4 of j2], the balls are assumed to be distributed along a hypersurface (for instance, 3d balls 
with centers periodically located on a plane). In the setting considered here, the rough idea is 
to extend the Stokes solution to the lower half-space by appropriate symmetry: our problem 
is then reduced to the homogenization of the Stokes equations across a planar network of 
patches. Hence, the ideas of [2], devoted to a planar network of balls, essentially apply. They 
are based on the construction of correctors and the method of oscillating test functions. We 
start with 

Lemma 6 (Existence of correctors). Assume that a e < e 2 . For every e > 0, there exist 
W e = (W&)l<iJ<3 e H 1 ^) 9 , <f = (Qj)i<j<3 G L 2 (nf, supported in T 2 x [-e/2,e/2], which 
satisfy the following properties: 

(i) W e — 1 weakly in H l (£l), q e — weakly in L 2 (f2); 

(ii) for every j = 1 . . . 3 ; ]T\ diW^ = in fi; 

(Hi) for 1 < i,j < 3, = Wf 3 = on T 2 x {0}, and for 1 < i,j < 2, Wfj = on 
T x {(.}.• 

(iv) For every <p G C°°(Q) 5 , every ip G H 1 ^) 3 and every sequence tp e G H 1 ^) 3 satisfying 
the boundary conditions 

ip e = on (T e x {0}) U (T 2 x {1}), ^3 = on T 2 x {0,1}, (3.3) 

and converging weakly to ip in H 1 ^) 3 , the following relation holds: if lim e _>o a t /e 2 = 
Cq > 0, then 

hm V ( f VWg • Vtffa - [ d^fq^A = -C [ Mo^ h ■ <p h . (3.4) 
i<S<3 Jn J Jt2 X {o} 



where Mo G 7^2(1^) is the symmetric definite positive matrix given by formula (3.19). 

Proof of Lemma\^ This lemma is the analogue of Proposition 4.1.6 in [3] (see also section 2.3 
in [2]). As mentioned before, we do not claim any major novelty in the proof. Nevertheless, 
with regards to quantitative aspects, notably the exact expression of the slip matrix Co Mo, 
we feel necessary to reproduce its main steps. 

The starting idea is to consider a base flow (W, q) in the vicinity of T, which, after proper 
rescaling, will describe accurately the corrector behavior near a single patch. We shall then 
truncate it and periodize so as to obtain an appropriate global corrector. Namely, we introduce 
the solution (W, q), with W = (Wij)i<i,j<'i-> Q = (%)i<i<3; °f the following problem: 

-AWij + diqj = in R 3 + , 1 < i,j < 3, (3.5) 

3 

diWij = in R%, 1 < j < 3, (3.6) 

i=l 

completed with the boundary conditions 

Wij = 8 i:j on T x {0}, 1 < i,j < 2, (3.7) 
W i3 = W 3j = on M 2 x {0}, 1 < i, j < 3, (3.8) 
d 3 Wij = on (K 2 \ T) x {0}, 1 < i,j < 2. (3.9) 
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as well as lim| a ,|_ 5>00 W = 0. Of course, for j = 3, we have Wis = and qs = 0. For j = 1, 2, 
the existence of a unique weak solution 

(W. jiqj ) G (Z) 1 ' 2 ^)) 3 x L 2 oc (M 3 )/M 

follows from Lax-Milgram theorem. Following |9j, we remind that L> 1,2 (IR^) is the closure of 
D(Rj) in fl" 1 ^). 

Asymptotic behaviour ofWij,qj. For j = 1...3, we extend qj, Wij and W2j into even 
functions of X3, and W$j into an odd function of X3. We obtain in this way solutions of 
the Stokes equations on M 3 \ T. Proceeding exactly as in j2j p. 255] , we obtain the following 
asymptotic expansions 



W.j (a?) 
qj(x) 



1 / F, 



8tt \\x 



+ 



(x ■ Fj)x 



1 + 1 



47T X 



+ 0( 1 — jp) as \x\ 



as |x| — >• 00, 



00. 



(3.10) 
(3.11) 



In formulas (3.10)-(3.11 ), the notation Fj corresponds to the drag force, which is defined by 
(here, n+ := e 3 , n_ := -e 3 ): 



Fa 



Tx{o+} 9n + 



-2 / 3 W.j, 

'Tx{0+} 



Tx{0"} 



<9n_ 



+ 



Tx{0+} 



Tx{0"} 



(3.12) 



Construction ofW e and q € . Using the extended W and q, defined in the whole of M 3 , we 
can then proceed exactly as in [2J, [3] to construct the correctors W e and q e . Therefore, we 



consider the following decomposition of P k (see definition (3.1)): 



Pj=CiUDlUKi, 

where C k is the ball of radius e/4 centered in the cube, D k is the ball of radius e/2, with same 
center, perforated by C k , and K k is the remaining part of the cube, that is K k = P k \ Dl (see 
Figure [3|). We denote by c| the center of cube P k . In each part of the cube, we define W.„- 



and qj as follows: 



") w r<e (V^-AW 6 , =0 . n 
1 , , / ' ~ c % \ V;egC ^> divW^io mD 



W e = 
7 _ ini^ 
d = fe 



Moreover, we impose J DE gj = and W% £ H 1 ^) 3 (so that there is no jump of W e across 
dDldCl). 

Estimates on W e and q e . We use again the decomposition P k = C k U (Pj: \ C|). The 
estimates in C k follow from the asymptotic expansions (3.10)-(3.11 ) and a scaling argument: 
for every e > 0, 

\\VW!j\\l HCV <Ca e , IkjH^ce) <Ca € , \\W% 2{c%) <Cale, 
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■< >► 

e 

Figure 3: Each cube P k is decomposed into a union of subdomains Cjr,, D e k and which 
separated by spheres of radius e/4 and e/2 centered in the cube. 



where C > is a constant. To treat the remaining part P| \ C|, we use a properly rescaled 
version of standard estimates for the homogeneous Stokes equations: basically, the L 2 norm, 
resp. H 1 norm of the solution is controlled by the I? norm, resp. H 1 / 2 norm of the boundary 
data (see for instance [E]). Since the velocity fields W e j satisfy the following pointwise 
asymptotics as e vanishes 

W e j = O on dC e k n 8D%, VW e j = O (J) on dC% D dD e k , 

using a scaling argument, we obtain the following estimates 

2 a 2 a 2 e 2 2 

ll VW -jll L 2(pe\qp ^ C -f> H^lli 2 (P|\Cf) - C ~^> \\ W ^L^(Pi\cl) - Ca ^ ( 3 - 13 ) 

for a given constant C > 0. Since < a e < e, we deduce 

HVW^-llia^.) < Ca e , \\q]\\^(P D < Ca e , < Ca 2 e. 

As a result, summing over k £ [|0, e _1 |] 2 , we obtain the following asymptotics as e vanishes 



|Vink (n) =0(5), lk e HL^)=0(5), ||^||f 2(n) = Off). (3.14) 



Conclusion of the proof. Let G C°°(0) 3 , i\> E F 1 ^) 3 and let ^ G i? 1 ^) 3 be a 
sequence of vector fields satisfying the boundary conditions (3.3), and converging weakly to 
ijj in H 1 ^) 3 . In the sub-critical case a e <C e 2 , the asymptotics (3.14) imply that 

W e -> strongly in ^(O) 9 , g e -> strongly in L 2 (S1). 

Consequently, the following relation holds: 

6 l<j,j<3 Jn ' 
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Thus, relation (3.4) holds with M = 0. 

In the critical case lim^o = Co > 0, we define Cl = T 2 x (— 1,1), and we extend (f>j and 
ijjj into even functions of x$ on for j = 1,2, and <j) 3 and "03 into odd functions of x 3 . First, 
asymptotics (3.14) imply that W e is bounded in H 1 , and therefore converges weakly in H , 
up to a subsequence. Since W e vanishes in L 2 (Q), we obtain VW e — 1 weakly in L 2 (J7) 9 . 
From (3.14), we also infer — ^ weakly in L 2 (£l), and thus the following identity holds for 
every 1 < i < 3, 1 < j < 2 

Moreover, 



VWfj ■ V(V#j) - q e jdi(il)!<t>j) + o(l), as e -> 0. 



= E / « + v ^ • v (^<) " ^(«) 

+E / e _vw-.-v(^)-^(^D + E / E e ^-vM-#i(M £ )> 

where C^' 1 * 1 = C| fll|. In all sums, k ranges over [| , e 1 1 ] 2 - Using the estimates (3.13), we 
infer that 



< c 



at 1 



e e z 



1/2 



< 1. 



At this stage the proof differs slightly from the one of [3], because of the mixed boundary 



conditions at x 3 = 0. Indeed, since (W t ,q t ) satisfies the Stokes system in C^, we have 



c, 



c; 



E,± 



<9n 



where n denotes the outer normal to the set C fc ' . In particular, due to the symmetry 
properties of W e , q e , (j), ip e , there holds 

f vwg ■ v(Md - [ e + fMMt) + I e _ ™h ■ v (Mt) - [ ei _ tfMMi) 



dCt 



cin{z=o} 



By definition of W, d 3 W^ = on (C e k n {x 3 = 0}) \ (T e x {0}). On the other hand, since 
■0 e satisfies (3.3), ipf = on T e x {0}. Therefore, the r.h.s. reduces to the integral on <9Cj;. 
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Using the asymptotic expansions (3.10)-(3.11 ) and the expression of Wfj,qj in we obtain, 



on 



dC 



dW. 



'j 



dn 



2 r, 6 r, 

—Fij + —Fj ■ n k nk ■ ei 

TT TT 



+ ( J) eRlj, 



where n k is the unitary normal vector to the ball C\ and R\^ is a function of x, satisfying 
RfJx) = O(l) as e — )• 0, uniformly in x and k. This leads to the following decomposition 



2 6 

*J "t J ' ^k^k ' &i 
TT 7T 



v 6 7 , J9C, e 



Let <5gc^ be the unit mass concentrated on dC^. Let us recall the following results, which are 
proved in [3] (see Lemma 4.2.1): the radius of the ball C| being equal to e/4, the following 
convergence hold^] 



TT 

16 



^T 2 x{0} strongly in H 



^n k n k ■ eiSgc* -> ^i^x{Q} strongly in H 1 (f2). 



(3.15) 



Let us write 



H-l(Q)xHl(a) 



(3.16) 



consequently, since — > 4>jipi weakly in H l (&) and % — > Co, we obtain 



2 6 

-^27 "I - 7 ' ^k^k ' &i 
TT TT 



T 2 x{0} 



(3.17) 



Moreover, since is uniformly bounded in L°°(J7), pondering on the comparison principle 



of [21 Lemma 2.3.8] we deduce from (3.15) that 

e Rtj5dC* k -> strongly in tf" 1 (O) . 



Using (3.16), we obtain the following convergence: 

eRfj^jipl — > as e ->■ 0. 



E 



k - dc t 



(3.18) 



1 Notice that in the paper of Allaire, the periodicity of the pattern is 2e, rather than e as in the present 
paper. Hence the constant in front of the Dirac mass in the right-hand side is tt/16, rather than 7r/64 for the 
first line, and 7r/48 rather than 7r/192 in the second line. 
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Gathering the convergence results (3.17)-(3.18), we obtain relation (3.4), where the matrix 
Mq is defined by 

' (3-19) 



Mi 



0,ij 



F ij} F.j given by (3.12) 



There only remains to prove that the matrix (Fij)i<i,j<2 is positive definite. To that end, 
we go back to system (3.5)-(3.9). We multiply by W.i the system satisfied by W.j, and we 
obtain 



F, 



'j 



-2 [ dsWij 

JTx{0} 



VWa ■ VW.. 



In particular, for all r] G 



ViVjFij = 2 I \V{r]iW.i + mw. 2 )\ 2 > o, 

l<i,j<2 

and the right-hand side above vanishes if and only if rjiW.i +112W.2 = a.e. in Ri. In view of 



6. 



the boundary conditions (3.7), this implies rji = 7/2 = 0. This concludes the proof of Lemma 

□ 
3 s 



To complete the proof of Theorem [l we rely on Lemma 6j as follows. Let <fi G C°°(iT 
satisfying the no-slip condition = on the upper boundary T^<{1}, and the non-penetration 
condition 4> 3 = on the lower boundary T 2 x {0}. Let W £ G H 1 ^) 9 , q e G L 2 (0) 3 be 
the sequences introduced in Lemma [6} We define the following test functions for the weak 
formulation associated to system (2.1 )- (j2.3h : 



<f = (I 3 - W*)<f>, r e = <f </>, 
where I3 is the identity matrix in A^3(M). We deduce the following relation: 



Vu € : V4> € - / p e dw 
n Jn 



n 



I 

Jn 



r e div u e = 0. 



(3.20) 
(3.21) 



Since W e converges weakly to in H (fi) , and strongly to in L (Q) , we readily obtain 

/ Vu e : V«f = / Vn e : V<f> - [ VufVW^ - / V<V<^- 
Jn Jn Jn Jn 

= Vu:V(j)- Vu\VW^(j)j + o{l), ase^O, 
Jn Jn 



— / p e div 
'n 



/ p e (I 3 - W € ) : V(f> = - / p div (j) + o(l), as e -»■ 0. 
Jn Jn 



Consequently, summing relations (3.20) and (3.21), we deduce the asymptotic relation 



/ Vu : V0 — p div 
In Jn 



/ pdiv<p- / VuiVW-j^j + / q € <p div u e = o(l), as e 0. 
Jn Jn Jn 
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Applying Lemma [6] with ip = u and ip e = u e , we obtain the following relation: 



Vu:V<f>- / p div 4> + C / M u -0 = 0, 



(3.22) 



where the matrix Mq € is defined by (3.19). Since relation (3.22) holds for every test 

function tfi, this proves that u = uc Q M - 



4 Asymptotic study of "riblet" designs 

This section is devoted to the proof of Theorem [2] In the case of riblets, we recall that T € is 
invariant by translation in x%. Since / = (/i,/2,/s) is also independent on the x\ variable, 



the solution (u e ,p e ) of system (2.1)-(2.2)-(2.3) depends only on (x2,xs). As a result, the first 
component of u e satisfies: 



d 3 u\ 



-A 2 , 3 «i 

0on(Tx{0})\(nT £ ), u\ 



h inTx(0,l), 
onTx{l}, 

o on nr e , 



(4.1) 



where V2,3 and A2,3 stand for the gradient (resp. the Laplacian) with respect to the (£2, £3) 
variables, T 1 = M/Z and where we have denoted II the projection operator defined by 
n(xi,X2,0) = (x2,0). In the same fashion, (u2,u^),p € satisfy the following Stokes problem: 



-A 



2.3 



■it 
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+ V 2 ,3P £ 



3, 
V 2 ,3- 



inTx (0,1), 



■(/■■: 



II. 



2.3 



a 3 ^ = 0on(Tx{0})\(nr e ), n 2 . 



x 3 =0 



inTx (0,1), 

onTx{l}, 

on T x {0}, 

on nr e . 



(4.2) 



Hence, the original 3d problem reduces to the study of two independent systems (with a 
Laplace and a Stokes equations), set in the 2d domain T x (0, 1). This change from a 3d to 
a 2d setting explains the change of scalings between Theorem [l] and Theorem [2] 



To handle the Stokes equations (4.2), we proceed like in the previous section: in short, 



we adapt the homogenization techniques of [3], dedicated to the Stokes flow across a 
periodic network of balls, set along an hypersurface. As mentioned before, the difference is 
the dimension of the domain. One must this time consider the 2d results of [3], about periodic 
network of disks along a line. For brevity, we do not give further details. We eventually obtain 
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the following limit system. 



-A 



2,3 



+ V 2 , 3 p = ( f A in T 1 x (0, 1) 
= in T 1 x (0, 1) 



u 2 ,3 = on T 1 x {1}, 
u 3 = on T 1 x {0}, 



(4.3) 



2tt 



d 3 u 2 = — on T 1 x {0}. 



As regards the Laplace equation (4.1), the idea is exactly the same. Actually, the situation 



is even simpler, and has been analysed for a longer time. Namely, one may start from the 
work of Cioranescu and Murat [7], instead of j2j section 4]. Again, we leave the details to the 
reader. In our setting, the limit system is 



-A 2) 3«i = /i in T 1 * (0,1), 
ui = onT 1 x{l}, 



7T 



(4.4) 



dzu! = ^—u x on T 1 x {0}. 



We deduce from systems (4.4) and (4.3) that u = UM riblets , M r n,i ets being given by (2.8). 
The sub-cases where / = e\ or / = e 2 follow easily. 

5 Numerical simulations 



This section is devoted to simulations of system (2.1)-(2.2)-(|2.3h. For simplicity, we shall 



restrict to constant source term (average pressure gradient), say 

/ = 2e, e e span(ei, e 2 ). 

The idea is to recover numerically the scalings for the slip length given in Theorems 1 and 
2. However, to observe significant slip implies to consider very small scales: patches of size 
less than e 2 , in a grid of side e. This forbids direct computations. To overcome this difficulty, 
we shall rely on a boundary layer approximation of the Stokes flow. Such approximation, 
often implicitly used in physics papers, has been fully justified in the context of wall laws: 
see [T21 IHJ HJ among many others. 

The starting point is to write the exact solution u e as 

u e (x) = u p (x) + ev e (x/e) 
where u p is the reference Poiseuille flow, satisfying (2.1) with Dirichlet condition at both 

u p [x) = -x 3 (x 3 - l)e. 



planes. Remind that 
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Hence, v e = v e (y) satisfies 

-Av + Vp = 0, in T 2 x (0,1), 
div v = 0, in T 2 x (0,1), 

v = 0, y 3 = e~\ (5.1) 
v 3 = 0, y 3 = 0, 
v h = 0, x{0}, d y3 v h = -e, y G e _1 (T e ) c X {0} 

Note that no approximation has been made so far. It is then tempting to put the roof y 3 = e _1 
at infinity replacing T 2 x (0,e _1 ) by T 2 x M_j_. However, it is well-known that the resulting 
problem is over deter mined. Namely, the boundary layer field v e ' 1 satisfying 



div v = 0, in fx 
^3 = 0, y 3 = 0, 



-Av + Vp = 0, in T 2 x R + , 

M+ ' (5.2) 
v h = 0, yee- x T e x{0}, d ys v h = -e, y € e~ l {T e ) c x {0} 
has constant horizontal average: 

e.oo / e.bl / \ 7 7 

v h '■= / v h vy)dy\dy2 
with respect to y 3 . More precisely, it can be shown that 



exponentially fast as y 3 goes to infinity. Furthermore, by linearity of (5.2), one may denote 
v h°° = y e '°° e for a two by two matrix V e '°°. Then, one can show that V e '°° is symmetric 
positive definite, with 

f \Vv e ' bl \ 2 . 

JT 2 xR + 



Note that everything depends on e, through the rescaled domain e T e in (5.2) 



To correct the "boundary layer constant" at infinity, one must add a macroscopic Couette 
flow. One ends up with 

u e ps u p (x) + ev e ' h \x/e) - ex 3 {V e '°°e, 0) 
Averaging in the small scale, we find 

u t,h \x 3 =o ~ eV e '°°e, d 3 u e \ x . i=0 ps d3^ P U 3 =o ps e. 
We end up with the approximate boundary condition 

u\ = eV e '°°d 3 u e h at x 3 = 0. (5.3) 

On the basis of the previous reasoning, one can implement the following strategy for the 
numerical computation of the slip length: 
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• Compute numerically (say with e = e\ and e = e^) the solution of (5.2), in order to 
determine the matrix V e,oc . 

• Check for the asymptotics of eV e,0 °, for various shapes and sizes of the no-slip zone 
T € . This allows to make the comparison with theoretical results of Theorems 1 and 2. 



Indeed, sending e to zero in (5.3) yields 



u h = lim ( e y e '°°) 8 3 u h at x 3 = 0, (5.4) 

e->0 

so that the matrix M in the theorems satisfies M~ 1 = lim (E _ 5 .o (eV e,co ). 
Numerical approximation of the matrix V e '°°. In the numerical simulations, we will 



solve the system (5.2) associated to different shapes of the no-slip zone T e : circular or rect- 
angular patches, and riblets parallel or orthogonal to the flow. Let us first notice that for 
such configurations, the matrix V €,co is diagonal. Indeed, since the domain e~^T e is sym- 



metric with respect to the axis {2/2 = 1/2}, if we denote by v the solution to system (5.2) 
with e = e\, then the vector field v* defined by v*(yi, y2, y 3 ) = Vi(yi, 1 — 2/2,2/3), for i = 1,3, 
and by 2/2, 2/3) = —"2(2/1, 1 — 2/2, 2/3), is also a solution. By uniqueness, we deduce that 

"2(2/1, 1 - 2/2,2/3) = -^2(2/1,2/2,2/3) for a.e. (2/1,2/2,2/3) G T 2 x R+, which yields V e >°°ei • e 2 = 0. 
By symmetry of V e '°° , we obtain also that V e '°°e2 ■ &\ = 

Consequently the boundary conditions satisfied by the horizontal components of the ap- 



proximate solution to system (2.1)-(2.2)-(2.3) on X3 = 0, simply writes: 



u\ = e(y e '°°ei • a) d 3 u e i at x 3 = 0, for i = 1, 2. (5.5) 

In the rest of this section, for i = 1,2, the quantity V e '°°ei ■ will be refered to as the average 
slip length associated to our problem, in the direction e%. 



To compute an approximate value of the average slip length associated to system (5.2), 
we consider a truncated domain T 2 x (0, H), for a given H > 0, and we introduce the solution 
w to the following problem: 

-Aw + Vq = 0, in T 2 x (0,#), 
div w = 0, in T 2 x (0,H), 
d y3 w - q e 3 = 0, y 3 = H, (5.6) 
^3 = 0, 2/3 = 0, 
w h = 0, i/er^xjO}, d y3 w h = -e, y G e _1 (T e ) c x {0} 



We solve problem (5.6) by a finite element method. We use P2 elements for the velocity and 
P\ elements for the pressure. The three-dimensional mesh of the fluid domain T 2 x (0, H) 
is obtained by a constrained Delaunay tetrahedralization. The incompressibility condition is 
treated by a Lagrange multiplier (see [10] , [11] ] 



Given two approximate solutions wl pp ,w^ pp of system (5.6), associated respectively to 



e = e\ and e = ei, we define the numerical approximation V a p P of the matrix V e '°°, by the 
following formula: 



v aw &i ' e i := J J2 w app(yi>yi> H ) • e 3 d Vidy2, for i,j = 1, 2. 
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By analogy with formula (5.5), for i = 1,2, the approximate average slip length in direction 
e{ is then defined by Vapp e-% • &%■ 

Finally, we introduce the solid fraction (jf s , which is defined by the relative area of the 
no-slip zone T e in the elementary square of size e (or equivalently, by the area of the rescaled 



no-slip domain e T e ). Using definitions (2.4)-( |2.5[ ), (jf s is given by the following expressions: 



in the case of patches, 
in the case of riblets, d 



\T\, where \T\ stands for the area of the domain T; 



\I\, where |/| stands for the length of the interval I. 



Notice that system (5.2) is completely determined by <fi e s and by the domain T (in the case of 



patches) or the union of intervals I (in the case of riblets). 



Computation of the average slip length, in the case of patches. In the case of 
patches, we have plotted Vapp £\ • ei against 1/ considering circular and squared patches 
(see Figure [4]). We observe that the dependency is affine, and a linear regression gives the 
relation Vapp^ei ■ e\ a/y^ + /3, with a = 0.322, /3 = —0.429 in the case of the disk, and 
a = 0.311, /3 = —0.422 in the case of the square. Note that these coefficients are very close 
to the ones obtained by Ybert et al. JIT]. Consequently, since lim e _>.o <t>% = 0, 



a 



as e 



0. 



(5.7) 



To compare this numerical result with the theoretical result given by Theorem [TJ let 
us consider the critical case a f /e 2 -> C > 0. In that case, there exists a two by two 
matrix Mo, depending on the pattern T, such that lim^o eV 6 ' 00 = ^M _1 . For circular 
or squared patterns centered in the unit square, as observed above, the matrices V e '°°, and 
consequently the matrix Mo, are diagonal. Moreover, since these patterns are invariant by 
a rotation of angle tt/2, one can easily see that the corresponding matrix V e '°° satisfies 



V €,co ei ■ e\ = V € '°°e2 ■ e2- Consequently, there exists Ao > such that Mo 
following relation holds: 

limey e '°°ei • ei ' 



Ao o 
Ao 



, and the 



Co Aq 



Besides, using the definition of d>\ in the case of patches, the asymptotic relation (5.7) yields 



J^meV^ei • e x 



o 



Thus, the numerical value of the slip length a/(Co\/\T\), that can be deduced from the 
asymptotic behavior (5.7) in the critical case, is consistent with Theorem [T] The coefficient 
of the matrix Mo can be approximated by Ao ~ \f\T\/a. 

We notice that the results concerning the sub-critical and super-critical cases can also 
be retrieved, at least formally, from relation (5.7). Indeed, since e/-^/0| = e 2 /(a, 
obtain in the sub-critical case 



lim f 



le->0 eVappei ■ ei 

infinite slip length in the e\ direction, that is, a perfect slip condition. In the same manner, 
in the super-critical case, we obtain lim e _s.o tVapp^ei ■ e\ = 0, which corresponds to adherence 
in the e\ direction. 



ey(a e ^\T\), we 
+oo, which corresponds formally to an 
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Computation of the average slip length, in the case of riblets. In that case, exact 
computations are available in the literature, that give the average slip lengths in the e\ and 
e2 direction as a function of the solid fraction <j> e s (see for instance [14]): 



ex ■ ex 



In 



7T 

cos l — (1 



e-i ■ e2 



In 



cos 



7T 



(l-#) /(2tt). (5.8) 



We have plotted in Figure [5] the computed value of the average slip lengths Vap^ex • ex and 
Vapp 3 ^ • e2, against as well as the exact values defined by formulas (5.8). We observe that 
the numerical values are close to the expected ones. 



Once again, formulas (5.8) and the numerical behavior of the average slip length shown in 



Figure [5| are consistent with the theoretical results of Theorem [2j Indeed, in the critical case 
lim e _;.o — eln(a e ) = Co > 0, using the expression <ff s = (a e \I\)/e, one obtains by a straightfor- 
ward computation that eln [cos (§(1 — 0s)) 1 ~~ ^ — C) as e — K Consequently, the slip length 
in the directions ex and e2 are respectively given by 



limey e '°°ei -ei 

e->0 



Co 



7T 



lim eV e '°°e2 • e 2 



Co 
2tt' 



Influence of the shape of the no-slip area: comparative results. In order to provide 
a comparison between the efficiency of patches and riblets in terms of slip length, we consider 
the slip length in the direction of the constant pressure gradient / = 2ej, with i = 1 or % = 2. 
For circular or squared patterns, the average slip length is given by V e '°°ex - ex- In the case 
of riblets, we consider two configurations of physical interest: 

• riblets parallel to the flow: / = 2ei, the average slip length is defined by V t '°°ex ■ ex', 

• riblets orthogonal to the flow: / = 2e2, the average slip length is V e '°°e2 ■ e 2 . 

The results are plotted in Figure |6] As stated in Remark [3j page[5j these numerical results 
confirm that the riblets parallel to the flow are not necessarily optimal. Indeed, if the solid 
fraction (jf s is small enough, say <f>\ < 0.1, the circular or squared patches produce a superior 
slip length. 

To estimate the influence of the shape of the pattern on the slip length, we have con- 
sidered families of rectangles of fixed area (jf s , that are centered in the unit square. For 
(f> e s = 0.01, 0.04, 0.09 we have computed the average slip length V^ex • ei, in the direction ex, 
associated to each of these rectangular patterns. The results are plotted in figure [7j against 
the dimension L of each rectangular pattern, in the e% direction. For each solid fraction <f>l, 
the extremal values associated to L = <jf s and L = 1, correspond respectively to a riblet 
orthogonal to the flow, and parallel to the flow. 

We notice that, for each family of rectangular patterns of fixed area, the riblet orthogonal 
to the flow provides always the smallest average slip length. As already mentionned, the riblet 
parallel to the flow is not optimal, especially for small values of the solid fraction <j>l = 0.01, 
(j)g = 0.04. In that cases, the curves present a unique maximum, and the associated optimal 
size L of the rectangle is slightly superior to the size v/0| of the square of same area. For these 
values of the solid fraction, the optimal rectangular pattern will present a certain anisotropy 
in the direction of the flow. 
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6 



Figure 4: Numerical value of the average slip length Vap^ei ■ e± plotted against for 
circular patches and squared patches. 




Figure 5: Numerical values of the average slip lengths Vap^ei ■ e\ and Vapp &2 • e2, plotted 
against <ff s , in the case of riblets. The dashed lines represent the exact value of the average 
slip lengths, defined by formulas (5.8). 
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5 • 10- 2 0.1 0.15 0.2 0.25 



Figure 6: Average slip length in the direction of the flow, in the case of circular patches, 
squared patches and riblets parallel and orthogonal to the flow. 




0.2 0.4 0.6 0.8 



Dimension of the rectangle along e\ 

Figure 7: Numerical values of the average slip length Vap£e\ ■ e±, produced by rectangular 
patterns of given area <j) e s G {0.01,0.04,0.09}, and plotted against the dimension L of the 
rectangle in the e\ direction. The extremal value L = (j) e s (resp. L = 1) corresponds to the 
riblet orthogonal to the flow (resp. parallel to the flow). 
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Appendix: proof of inequality (3.2) 



To obtain inequality ( |3.2[ ), it is enough to prove that for every k G [|0, e 1 |] 2 , 

\u e \ 2 < 77(e) / |Vu e | 2 VO < e < e . (5.9) 



'Sjrx{0} J B 

A summation over k G [ 1 , e 1 1 ] 2 then leads to inequality (3.2) 



Let k G [|0, e 1 |] 2 . By rescaling the trace inequality in the half cube [0,1] 2 x [0, k], we 
obtain the existence of a constant C > such that 

|«f<C(e / |W| 2 + - I K| 2 V (5.10) 



5fx{0} \ ^ + 



To estimate the L 2 norm of u e by the 1? norm of its gradient, we adapt [H Lemma 3.4.1] 
to our bidimensional framework. We denote by B k the ball circumscribing the cube P k . Of 
course, the upper half-cube Pf' is contained in the upper half-ball Bf' . Moreover, since 
the model no-slip zone T contains a disk of radius a centered at the origin, each elementary 
no-slip pattern ek + T e contains a disk of radius a e a, centered in the square S k . Let B k be the 
ball of same center and radius, and B k + be the corresponding half ball. With this notation, 
we can write 

| u f < / | n f + / 1^12 



To estimate the contribution of the exterior part B k + \ B k + , we use spherical coordinates 
(p,4>,6) centered at point ek + (| , |, 0). The radius of B e k being equal to integrating 
along rays, we get for every r', r such that < r' < a e a < r < 



u\r^,9) = u e (r',cj>,9) + / d p u £ (p, <t>, 0)dp, 

Jr' 

which yields 

\u\r, 0, 9)\ 2 < 2\u e (r', 0, 9)\ 2 + 2 (£ d p u*(P, <t>, 0)dp\ . 

Multiplying last inequality by r 2 (r') 2 sin 8 and integrating on r' G (0,a e a), r G (a e a, ^^), 
(j) G (0, 27r), 9 G (0,7r/2), we obtain the inequality 

I e <2J e + 2K e (5.11) 

where the integrals I € , J e , K e are respectively defined by 



r'=0 J r=a fL a 



u(r,cf),9)\ r (r) sin# d9 dcj) dr dr' , 



J e = / / / / \u(r l ,(j),9)\ 2 r 2 (r') 2 sm9 d9 d<j> dr dr', 
K e = I ' I 2 [ [([ d p u e (p,<p,9)dp) r 2 (r') 2 sm9 d9 d(p dr dr'. 

Jr'=0 Jr=a e a J 9 J <h \Jr> J 
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By Fubini theorem, 

° (r'fdA ( f 2 
and by an analogous computation, 

J e = 

By Schwarz inequality, 



u r, 



I V sine d9 deb dr 



a 3 a 3 



.e\2 



U 



,3^3 a 3 a 3 



■a 



c|2 



B, 



u: 



d p u e dp I < 



1 

r> P 1 



-Ap 



p l \d p tf\ l dp ) < _ / p z \d p u e \ 2 dp, 



which yields 



K £ < 



r 2 dr 



r dr' 



Vu £ \ 2 ) < ~-^aia*e 
] 16 



3^/3 



B. 



iwi 2 . 



B, 



Consequently, inequality (|5.11|) leads to 

/ 

JbI:+\bI 



2 9^3 e : 

u < 



a e a I a 2 a 2 Js e,+ 



K| 2 + 



c]2 



*2'" 



(5.12) 



since u vanishes on n (R 2 x {0}), using Poincare inequality in a cylinder of height a e a, 
we obtain the following estimate: 



I K| 2 < a 2 a 2 [ 
JBV+ JbI 



c 1 2 



Injecting this inequality into estimate (5.12), we obtain: 

27\/3e 3 



\u £ \ 2 < 



i, \ B i 



a e a J B 



c 1 2 



and summing these two inequalities, we get 



I e 1 2 ^ I 2 2, 



27\/3e 3 
a £ a / J B 



\Vu £ \ 2 . 



Finally, using inequality (5.10), we obtain estimate (5.9), where 77(e) is defined by 

,s „( a 2 a 2 27>/3e 2 \ 

r/e) = C e+ ^ + , 

\ e a e a / 



and converges to as e — > 0, since a e < e and a e 3> e 2 . 



□ 
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